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We show that the a\pir Lagrangian is a decisive element for obtaining a good phenomenologi- 
cal description of the three-pion decays of the r lepton. We choose it in a two-component form 
with a flexible mixing parameter sin (9. In addition to the dominant a\ — > Tip intermediate states, 
the ai — > na ones are included. When fitting the three-pion mass spectra, three data sets are 
explored: (1) ALEPH 2005 u'ti + u~ data, (2) ALEPH 2005 tt'tt'V data, and (3) previous two 
sets combined and supplemented with the ARGUS 1993, OPAL 1997, and CLEO 2000 data. The 
corresponding confidence levels are (1) 28.3%, (2) 100%, and (3) 7.7%. After the inclusion of the 
ai(1640) resonance, the agreement of the model with data greatly improves and the confidence level 
reaches 100% for each of the three data sets. From the fit to all five experiments [data set (3)], the 
following parameters of the ai (1260) are obtained: m ai = (1233 ± 18) MeV, F ai = (431 ± 20) MeV. 
The optimal value of the Lagrangian mixing parameter sin (9 = 0.459 ± 0.004 agrees with the value 
obtained recently from the e + e~ annihilation into four pions. 

PACS numbers: 13.35. Dx,13.25.-k,14.40. Be 



I. INTRODUCTION 

The dominance of the ai(1260) meson, hereafter re- 
ferred as oi, and its irp decay mode in the three-pion 
decays of the r-lepton is firmly established experimen- 
tally [1-16]. A convincing demonstration was provided 
by the ARGUS Collaboration [7, 16], who compared the 
distribution of unlike- and like-sign two-pion masses. The 
■wp intermediate state is the core of several models of the 
three-pion decays of the r-lepton [12, 17-26]. It should 
be mentioned that the a\ dominance in heavy lepton de- 
cays was proposed in 1971 [27], about five years before 
the r-lepton was actually discovered. 

The ai(1260) resonance, discovered almost fifty years 
ago [28] , plays an important role in many phenomena of 
the nuclear and particle physics. Its properties have been 
studied in many processes, but even its basic parameters 
are not very well known. The values of the ai(1260) res- 
onance mass determined from different processes or by 
different experimental groups often contradict one an- 
other. The same applies, even to a larger extent, to the 
ai width. Very little improvement has been achieved over 
the last thirty years, see Table I. The origin of those 
problems lies in the very nature of the a\ resonance with 
its short lifetime and large width. The usual definition 
of the mass and usual procedures for its measurement 
are not applicable. The ai mass and width enter the 
formulas for experimentally accessible quantities via the 
assumed form of the resonance propagator, which gen- 
erates a specific Breit-Wigner formula. Those formulas 
are further modified in different ways when modeling the 
dynamics of the processes in which the a\ participates. 
As a result, we do not have a unique definition of the a\ 
mass m ai and width r ai . In fact, every formula repre- 
sents a specific definition of m ai and T ai . Given this, it is 
not surprising that different models yielded different re- 
sults even when being applied to the same data. It would 



TABLE I. Thirty- year history of the basic ai(1260) param- 
eters as listed in the Particle Data Group Tables. Only the 
last edition and those in which a change appeared are shown. 



PDG 


m 01 (MeV) 


r ai (MeV) 


ai — > pit 


1978 [29] 


« 1100 


« 300 


w 100% 


1980 [30] 


1100 to 1300 


« 300 


dominant 


1982 [31] 


1275±30 


315±45 


dominant 


1986 [32] 


1275±28 


316±45 


dominant 


1988 [33] 


1260±30 


300 to 600 


dominant 


1990 [34] 


1260±30 


350 to 500 


dominant 


1992 [35] 


1260±30 


« 400 


dominant 


1994 [36] 


1230±40 


« 400 


dominant 


1998 [37] 


1230±40 


250 to 600 


dominant 


2000 [38] 


1230±40 


250 to 600 


seen 


2008 [39] 


1230±40 


250 to 600 


seen 



be natural to accept as the a\ canonical parameters the 
results of a model that best describes a broad class of 
data on various processes and from various experiments. 
Unfortunately, we are not in such a situation yet. 

The situation of the heavier meson states with J PC — 
1 ++ is a little unclear and none of them has found its 
place in the Summary Table of the Review of Particle 
Properties [39]. The first indication of the state with 
mass of 1.65 GeV and width of 0.4 GeV appeared al- 
ready in 1978 [40]. The later experimental evidence, 
which comes mainly from hadronic reactions, is summa- 
rized in [39], where this resonance is listed as ai(1640) 
and assigned the mass of (1647 ± 22) MeV and the width 
of (254 ± 27) MeV. The three-pion decay of the r-lepton 
is less convenient for studying the ai(1640) resonance 
(often denoted as a[ in what follows) because of fun- 
damental limitations due to the r mass which is not big 
enough to provide sufficient phase space for three-pion fi- 
nal states with the needed invariant mass. Nevertheless, 
the DELPHI Collaboration [11] performed the analysis 
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of the Dalitz plots for different 3-pion mass ranges and 
observed an enhancement that "could be explained by a 
decay mode of the r to a resonance of mass similar to or 
greater than the r mass which then decays to three pions 
through the intermediate state of a pion plus a particle 
of mass 1.25 GeV or greater." They interpreted this as 
an evidence for the a[. In our opinion, this observation 
need not signify the existence of the a[ . It may also be a 
decay of the ai(1260), produced with a larger than nomi- 
nal mass, into ir and p(1450). A more convincing proof of 
the a[ in the decay of the r lepton comes from the CLEO 
Collaboration [12]. They showed that adding the a[ term 
into the Breit-Wigner function improved significantly the 
agreement with the data. 

On the theoretical side, a radial excitation of the 
quark- antiquark system with a mass of 1.82 GeV ap- 
peared in a relativized quark model with chromodynam- 
ics of Godfrey and Isgur [41]. Its decay width into the np 
channel was calculated in the flux-tube-breaking model 
by Kokoski and Isgur [42] with result < 70 MeV (our es- 
timate is based on their Table II). The seminal analysis 
of Barnes, Close, Page, and Swanson [43] has shown that 
the experimentally observed dominance of the D-wave 
over S-wave [44, 45] excludes the hybrid meson nature 
of the a[ and confirms it as a radial excitation of the 
quark-antiquark system. 

A few I G (J PC ) = 1~(1 ++ ) meson states above the 
ai(1640) have been observed by a single group, mainly 
in the pp annihilation. They still need confirmation. For 
details, see [39]. 

Another important ingredient that defines a particu- 
lar model of the three-pion decay of the T-lepton is, be- 
sides the a\ propagator, the aipir vertex. The models as- 
sembled by different authors are based on different aipn 
vertexes. These are sometimes simply constructed as al- 
lowed combinations of the metric tensor and participat- 
ing four-momenta. A more rigorous way lies in deriving 
them from the interaction Lagrangians among the ax- 
ial, vector and pseudoscalar fields. Unfortunately, here 
the situation is unclear yet. Various theoretical concepts 
provide different effective Lagrangians [46] . This is prob- 
ably the reason why the model builders preferred trivial 
Lagrangians or ad hoc vertexes. However, recent articles 
[24-26] are different. Dumm, Pich, and Portoles [24] got 
their Lagrangian from the resonance chiral theory. Their 
work was revised in the light of later developments in 
[25]. Achasov and Kozhevnikov [26] used the General- 
ized Hidden Local Symmetry model. 

Several models of the three-pion decay of the tau lep- 
ton have been proposed. With some simplification one 
can say that each of them gives compatible results when 
applied to different sets of data, but the results of differ- 
ent models are incompatible. Also the agreement of many 
models with data (often verbally claimed as satisfactory) 
is poor when judged by usual statistical criteria. The 
most popular models were those of Isgur, Morningstar, 
and Reader (IMR) [19] and of Kiihn and Santamaria (KS) 
[20] . Other models were much less successful in fitting the 



data. As an example we recall the results from [7] , where 
the ARGUS Collaboration compared various models with 
their data. Using the % 2 's and the numbers of degrees 
of freedom (NDF) from their Table 4, we are getting the 
confidence level (C.L.) of w 10~ 4 for Bowler's model [47] 
and 2.2% for the model of Ivanov, Osipov, and Volkov 
[48]. The KS and IMR models look better with C.L. 
10.7% and 79.0%, respectively. However, in a later arti- 
cle [49] the ARGUS Collaboration used an enlarged set of 
data (integrated luminosity of 445 pb _1 against 264 pb -1 
in [7]) and found that the KS model is rejected on a 7.4 a 
level. The IMR model with parameters as given in [19] 
was incompatible with the data on the same level [49] . 

Up to now, the best results have been obtained by 
the CLEO model [12] and by the model of Achasov and 
Kozhevnikov [26]. The former obtained, when fitting 
the CLEO 7r-7r°7r° data [12], C.L. of 54.6% without the 
a\ resonance and 88.2% with it. The latter fitted the 
ALEPH 7r~7r + 7r~ data [15] assuming two heavier axial 
mesons a[ and a'{ and got x 2 /NDF=79/102, which cor- 
responds to C.L. of 95.6%. Unfortunately, each of those 
two successful models has been applied only to one data 
set. 

The finding of an aipir Lagrangian that leads to a sat- 
isfactory description of the three-pion production in the 
tau decays would have important consequences for other 
areas of the high energy and nuclear physics. For exam- 
ple, the a\ resonance and its coupling to the pn system 
play important role in the evaluation of the dilepton and 
photon production rates from a hadronic fireball presum- 
ably created in the relativistic heavy ion collisions. The 
calculations performed so far, see, e.g., Refs. [50], have 
shown that the yield of electromagnetic signals strongly 
depends on the choice of the a\pn Lagrangian. Fixing 
its correct form is thus important for distinguishing the 
electromagnetic radiation of the Quark-Gluon Plasma 
(QGP) from the hadronic sources. 

The outline for this paper is as follows. In Sec. II we de- 
scribe our model and mention briefly its similarities and 
differences with other models. The experimental data 
used for testing our model and fixing its parameters are 
listed in Sec. III. Some details about our calculations 
and the results are presented in Sec. IV. We summa- 
rize our results and conclude in Sec. V. Two Appendixes 
contain technical details. The present work supersedes 
an earlier paper [51]. 



II. MODEL OF THE THREE-PION DECAYS OF 
THE TAU LEPTON 

In this section we present our model, which will be 
used for fitting the three-pion mass spectra of the de- 
cays t~ — > j/ T 7r~7r + 7r~ and r~ — > v T ir~ir°ir° . The basic 
information can be obtained by inspecting Figs. 1 and 
2. In addition to the standard a\ — > up intermediate 
states we include also the states in which the a\ couples 
to a pion and an / (600) (hereafter called a). The na 
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FIG. 1. Feynman diagrams of the r~ — > !/ T 7r~7r + 7r~ decay. 




FIG. 2. Feynman diagrams of the r — > ^ T 7r 7r°7r° decay. 

intermediate states improve the behavior of the differen- 
tial decay width at small masses of the three-pion sys- 
tem and bring the difference between t~ — > v T ir~ir + ir~ 
and t~ — > ^ T 7r~7r 7r° decays. Their importance has been 
pointed out by the CLEO Collaboration [12]. 

To show what is specific for our model, what differs it 
from other existing models of the three-pion decays of the 
tau lepton, we have to provide more information. This is 
done in the following subsections. 



example, [8]. 

To maintain both the flexibility and the correspon- 
dence with the effective field theory, we use a two- 
component Lagrangian of the aipir interaction in the 
form 

C aifm = (Ci cos 6 + C-2 sm6), (1) 

where 

& = A" • (V„„ x d"P) , 
C 2 = V„„ • (0"A" x P) , 

and V M „ = <9 M Vj, — c^V^. The isovectors A M , V M , and 
P denote the operators of the a\, p and 7r fields, respec- 
tively. 

Our Lagrangian differs from that derived by Wess and 
Zumino, see Eq. (67) in [52], only by notation. We will 
consider the mixing angle a free parameter that has 
to be determined by fitting the experimental three-pion 
mass distribution. For each 8, the coupling constant 
9a t piv can be determined from the a\ — > pir decay width. 
The Lagrangian (1) implies the following aipir vertex 

xa » = {cos9 [p> , M ga ^ 

where p's denote the four-momenta of the corresponding 
mesons (incoming a\, outgoing p and ir). 

Lagrangian (1) has recently been used [53, 54] in a 
model of the electron-positron annihilation into four pi- 
ons. Value of the mixing parameter sin 9 was obtained 
by fitting the excitation function (dependence of the an- 
nihilation cross section on the invariant collision energy). 
From the 7r + 7r~7r + 7r~ channel the value of 0.460 ± 0.003 
has been obtained [53]. In [54], a combined fit to 
7r + 7r~7r + 7r~ and ir + ir~ir°Tr° channels has provided the 
value of 0.466 ±0.005. 



A. Phenomenological a\p-K Lagrangian 

The interaction Lagrangian among the oi, p, and pion 
fields implies the form of the aipn vertex in the Feynman 
diagrams. But sometimes a vertex is postulated that can 
hardly be related to any effective Lagrangian. In the 
literature, one can find several prescriptions for the aipir 
vertex used in the calculation of the decay rate of the tau 
lepton into three pions and neutrino. The simplest one is 
X a ^ oc g afl , where index a (p) couples to the a\ (p) line. 
It can be derived from the interaction Lagrangian among 
the ai, p, and n fields without derivatives. It was used, 
e.g., in Ref. [17]. On the opposite pole of complexity is a 
two-component vertex used in the IMR model [19]. Both 
its components are transversal both to the ai and p four- 
momenta. The relative weight of the two components can 
vary, what gives the IMR model more flexibility. This is 
probably the main reason why this model sometimes fits 
the data a little better than the KS model [20], see, for 



B. Other effective Lagrangians and their 
parameters 

The Lagrangian describing the interaction of the a\ 
triplet with strange pseudoscalar and vector mesons is 
not required when calculating the amplitudes of the 
three-pion decays of the r, see Figs. 1 and 2. It is needed 
for evaluating the strange channel contribution to the to- 
tal decay width of the ai(1260). The latter enters the a\ 
propagator discussed below. The Lagrangian is chosen 
in a form analogous to Eq. (1) 

C ai K*K = gai ^ K cos 6 + A sin 0) . 

with 

£[ = d u K^A^Kl v + H.c. , 
C 2 = K^d^A v K*^ + H.c. 
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Matrix notation is now used, in which 



C. Propagator of the ai resonance 



K 



K+ 
K° 



K* 



K*+ 
n 

K*° 



VV2( 



(a?)" 



and if*,, = d^K* 



—d^K*. As usual [55] , a particle symbol 
denotes the field operator which annihilates that particle 
and creates its antiparticlc. In the spirit of the SU(3) 
symmetry, we assume the same mixing angle 9 as in the 
aipn case (1). The coupling constant g ai K»K cannot 
be reliably extracted from the experimental data yet be- 
cause of conflicting information about the a\ — > KK-k 
branching fractions. 1 . We will therefore use the SU(3) 
symmetry relation 



9 ai K*K 



1 



2 



(2) 



In order to evaluate the amplitudes of the Feynman 
diagrams depicted in Figs 1 and 2, we also need to specify 
the interaction Lagrangian among the a\, n, and a fields. 
We write it in the form 



r 



gi (A" ■ c^P) S + g 2 (A" • P) d^S , 



where S is the operator of the a field. The Lorentz con- 
dition for the ai field implies that the amplitude of the 
decay a\ — > a + n is proportional to the difference 



SaiCTTT — .9i — 92- 



(3) 



In the r decay diagrams, where the off-mass-shell a\ reso- 
nance is represented by its propagator (5), also the terms 
proportional to 



h ai cTTT — .91 + .92 



(4) 



contribute. There is no way of inferring the g ai cnr and 
h ai cm from the hadron decay data. We will return to this 
problem later in this article. 

The interaction Lagrangian between the a and 7r fields 
is given by 

The coupling constant g a ^ could be estimated from the 
data on the a mass and width [56, 57]. But because this 
constant enters the amplitudes of the three-pion decays 
of the taon multiplied by g aicrT! or h aia - Kl which are both 
unknown, it does not have much sense. 



We choose an analytically correct form [18, 19] of the 
a\ propagator featuring the running mass M(s) and the 
energy-dependent total width T ai (s) in the denominator 



f v +p»p v /m 2 ai 



Ml(s)+im ai T ai (s) ' 



The following conditions should hold 



<«) = 

dM t ( 2 , 

r ai (ml 1 ) = r ai , 



0. 



(5) 

(6) 

(7) 

(8) 



where m ai and T ai arc the nominal mass and width of 
the oi(1260) resonance, respectively. The denominator 
in (5) is the boundary value of a function analytic in the 
complex s-planc (s = p 2 ) with a cut running along the 
real axis from the three-pion threshold to infinity. The 
running mass squared can therefore be obtained from a 
once-subtracted dispersion relation 2 with T ai (s) as input 



M a 2 1 ( S )-M 2 i (0)--p/ 



m ai T ai (s') 

9ml S'(S'-S) 



da' . 



(9) 



Symbol P denotes the Cauchy principal value. We have 
chosen the subtraction point at s = 0, instead of s = m„ 
as in [19]. The advantage is that the integrand in (9) con- 
tains just one singular point instead of two, what makes 
the evaluation more stable and much faster. The disad- 
vantage is that the condition (6) is not satisfied auto- 
matically and M 2 x (0) must be recalculated if m ai or any 
parameter inside T ai (s) changes. 

The following decays are considered when calculating 

r 01 («): 



ai — >• p + 7r — >• 3tt 

ai -> K*K, K*K -> KKtt, 

ai — > a + tt — > 37T, 



(10) 

(11) 
(12) 



where the mass of the decaying a\ is taken to be y/s. 
We neglect the interference between the amplitudes of 
(10) and (12) despite the identical final states. We argue 
that the decay a\ — > p-n proceeds in the S and D orbital 
momentum states, whereas a\ — > an in the P state. This 
argument is not entirely watertight because neither (10) 
nor (12) satisfies the conditions for being factorized as 
a two-step process [58]. Channel (11) is described by 
four Feynman diagrams, two of them have identical final 
states (e.g., K-K°ir° in the case of a 1 .) We checked 
that the interference can be safely neglected in this case. 

The Lagrangian between the vector (K*) and pseu- 
doscalar (K, ir) fields is chosen in a standard form 



1 See [39] and the discussion on p. 253 in [15]. 



2 In Refs. [12, 18] an unsubtracted dispersion relation was used. 
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with coupling constant Qk*Kti- Moreover, the empiri- 
cal widths [39] of the p(770) and if* (892) provide the 
ratio 



9 2 k 



K- 



Q 2 



0.883 ± 0.035, 



(13) 



where the error has been enlarged to absorb the difference 
between the charged and neutral if* (892). Value (13) is 
a little higher than the SU(3) value of 3/4. Relations 
(2) and (13) enable us to express the product of coupling 
constants squared acting in (11) multiple of 



G 



2 2 



(14) 



which determines the partial decay width of (10). We 
introduce the ratio 



x = 



G 2 



(15) 



the value of which is given by multiplying (2) by (13). 

As we have already mentioned, there is no way of get- 
ting the product g ai o-n9o-K-K from the data, as the partial 
decay width of (12) is unknown. We therefore proceed in 
another way. We define the parameter 



y 



G 



(16) 



If we insert the parameters x and y into the formula 
for the a\ total decay width, it becomes proportional to 
G 2 . So does the derivative of the running mass squared 
(9). When the condition (7) is applied, G 2 can be can- 
celed. With known x, the condition (7) thus becomes an 
equation for the unknown y 2 . As we neglect the possible 
interference between the Feynman diagrams containing 
the p with those containing the a, the sign of y is not 
essential and we choose y > 0. The dimension of y is 
(energy) 2 because the two Lagrangians that describe the 
decay a\ — >■ 3ir via pir have together three derivatives, 
while those via an just one. 

An important note concerns the hadron vertexes. The 
effective Lagrangian approach takes hadrons as elemen- 
tary quanta of the corresponding fields, ignoring thus 
their internal structure. As a consequence, the inter- 
action strength is overestimated at higher momentum 
transfers. To describe the interaction among participat- 
ing mesons more realistically, we explore the chromoelec- 
tric flux-tube breaking model of Kokoski and Isgur [42] , 
as it was done already in the IMR model [19]. Each 
strong interaction vertex is modified by the factor 



F(q) = cxp 



Up 2 



(17) 



where q is the three-momentum magnitude of a daughter 
meson in the rest frame of the parent one (virtual masses 
are taken in the intermediate states). In the original 
paper [42], the value (3 = 0.4 GeV/c was established. We 
will use this value in all our calculations, as we found 



that moving from it did not bring statistically significant 
improvement of the agreement of our model with data. 
A cutoff similar to (17) was used in the model by CLEO 
Collaboration [12]. Their parameter R = 1.2 corresponds 
to £ = 0.340. 

During the course of development of our model we tried 
various versions of the a\ propagators, from the most 
primitive one with the constant mass and width to the 
most sophisticated and best physically justified one (5). 
The best fit to data has been provided by the latter. 

When investigating the presence of the suspected ra- 
dial recurrence of the ai(1260), denoted as a[, we sup- 
plement the a\ propagator (5) with the term 



' m \ + im a'J 'a[(s) 



(18) 



where a is a complex parameter. We assume that the 
energy dependent total decay width r a ' (s) exhibits the 
same energy behavior as that of ai(1260) and write 



r 01 («) 1 



F<»i ( m a; ) i s tne assumed width of the a[ 



where r o / 
resonance. 

What concerns the relation to the previous models, our 
di propagator is closest to that used by the CLEO Col- 
laboration [12]. If we ignored the momentum dependent 
terms in the numerators of (5) and (18), we would recover 
their Breit-Wigner function. 



D. Propagators of the p, K* , and a resonances 

In order to calculate the amplitudes of the taon's three- 
pion decays we need also the p and a propagators. They 
play a role also in decays (10) and (12). In addition, the 
evaluation of the decay width (11) requires the knowledge 
of the K* propagator. 

We choose the propagator of both the charged and 
neutral rho resonances in the form 



-iG^(p) 



-g» v + p^p u /m 2 p 
M 2 {s)+im p Y p (s) 



(19) 



which uses the running mass squared M 2 (s) and the en- 
ergy dependent total width T p (s) from Ref. [59]. The 
denominator of propagator (19) is an analytic function 
in the s-plane with a cut running from 4m 2 . to infin- 
ity, as required by general principles. The real function 
M 2 (s) is calculated from T p (s) using a once-subtracted 
dispersion relation, which guarantees that the condition 
M 2 {m 2 ) = m 2 is satisfied. The condition 



dM p ( 2, 







is not fulfilled automatically and serves as a check that all 
important contributions to the total p-meson width T p (s) 
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have properly been taken into account. They include, in 
addition to the basic two-pion decay channel, the W7r°, 
K + K~, K°K°, and r]ir + ir~, which get open as the p 
resonance goes above its nominal mass. The structure of 
the participating mesons is taken into account by means 
of the Kokoski-Isgur form factor (17). 

The running mass description of the p propagator [59] 
differs from other approaches that appeared in the liter- 
ature [60-62] . Gounaris and Sakurai [60] considered only 
the two-pion contribution to the total width of the p° res- 
onance and ignored structure effects. The result is a sim- 
ple analytic formula, the main reason why their approach 
is so popular. Vaughn and Wali [61] took into account 
the strong form factor, but again ignored higher decay 
channels. Mclikhov, Nachtmann, Nikonov, and Paulus 
[62] included the K+K~ and K°K° channels, but did 
not consider the strong form factors. The running mass 
formalism [59] takes into account both the higher decay 
channels and the structure effects. 

The propagator of the if* (892) resonance is required 
only for the calculation of the decay rate (11), which 
contributes to the total decay width of the ai(1260) res- 
onance. It does not act in the three-pion decay of the r 
lcpton. It is chosen in a simpler form, with the constant 
mass and energy dependent decay width 



-g^ + p^p v /m 2 K , 
s - m\, + im K * F K , (s) 



(20) 



The decay width includes only the contribution from the 
K* — > K + 7r channel and is normalized to the nominal 



width Tk* at s 



The corresponding formula, 



taking into account also the Kokoski-Isgur form factor 
(17), is 



r*.(*) = 



in 



K ' 



Q(s) 



1 3 



q(m 2 K ,) 



F(q(m 2 K ,)) 



K* 



where q(s) is the momentum of a daughter particle in 
the rest frame of the parent K* with the mass y/s. The 
if* (892) is a narrow resonance and we experienced nu- 
merical instabilities when calculating integrals containing 
the square of (20). To get rid of problems, we have used 
the procedure described in Appendix B. 

Also for the a propagator we use the form with fixed 
mass and energy dependent width 



-iG a (p) = 



1 



s — ml + im l7 T <7 (s) 



where V cr (s) includes only the contribution from the two- 
pion decay channel and is equal to 



m„ 



s - Ami F {l{ s )) 
sum 2 , - Ami F (q( m D) 



As the current Review of Particle Physics [39] is not 
very specific about the /o(600) mass and width, we rely 



on the mutually compatible values obtained by the Fer- 
milab E791 Collaboration [56] and the CLEO Collab- 
oration [57], who both analyzed the D mesons decays. 
The results (in MeV) of E791 are m a = 478j^ ± 17, 
T a = 324 ^ ±21, whereas those of CLEO are m a = 
513 ± 32, r CT = 335 ± 67. We adopt the weighted aver- 
ages m a = 500 MeV and r CT = 329 MeV. 



III. EXPERIMENTAL DATA 



We will compare the calculated three-pion mass dis- 
tribution in the r~ — > ^ r 7r _ 7T + 7r _ and t~ — > z/ T 7r~7r°7r° 
decays with the outcome of the five experiments. 

(1) The ARGUS Collaboration [7] used the ARGUS 
detector at the DORIS II e+e~ storage ring at DESY 
and studied the r~ — > ^ r 7r _ 7r + 7r _ decay. Their back- 
ground and acceptance corrected three pion mass distri- 
bution is given in twenty-eight bins within the mass range 
0.425-1.775 GeV. 

(2) The OPAL Collaboration published their results on 
the three-pion-mass squared distribution in the charged- 
pion channel in two papers. The first of them [8] was 
based on the data collected with the OPAL detector at 
the CERN Large Electron-Positron Collider (LEP) dur- 
ing 1992 and 1993. We used it in our recent publication 
[51]. In this work we explore the updated version [9], in 
which also the data of 1994 were included. The three- 
pion-mass-squared plot is corrected for background and 
efficiency and consists of twenty-three bins with much 
smaller statistical errors than in [8]. 

(3) The r-lepton decay into three pions and neutrino 
was also investigated by the CLEO Collaboration at the 
Cornell Electron Storage Ring (CESR). Their results on 
the all-charged-pions channel still exist only in a prelim- 
inary form [14]. We can therefore use only the data on 
the 7r~7r 7r° channel [12]. The background-subtracted, 
efficiency corrected three- 7r mass spectrum is given in 47 
bins. 

(4,5) The ALEPH Collaboration at CERN LEP have 
published an article summarizing their results about the 
branching ratios and spectral functions of the r decays 
[15]. It is based on the data collected with the ALEPH 
detector during 1991-1995 but processed by an improved 
method. We use the tables of the corrected three-pion 
mass squared spectra both in the t~ — > ^ r 7r _ 7r + 7r _ and 
t~ — > ^ r 7r _ 7r 7r° decays, which are publicly accessible at 
the website [63] . The all-charged pion spectrum contains 
116 bins 0.025 GeV 2 wide starting at 0.225 GeV 2 . In the 
two-neutral pion case the spectrum starts at 0.2 GeV 2 , 
but we discard the bin centered at 0.2375 GeV 2 with a 
zero value and, comparing to its neighbors, an unrealisti- 
cally small error. We are thus left again with 116 bins. In 
both cases we ignore the correlation matrices among the 
errors in different bins and add statistical and systematic 
errors linearly. 
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IV. CALCULATIONS AND RESULTS 

To be sure that our results are free of programming 
errors, we have written two independent computer codes, 
one in C++ (M.V.), another in Fortran 95 (P.L.) and 
debugged them until they produced identical results. 

Wc found that the parity-violating term in the r de- 
cay amplitude influences the three-pion-mass distribution 
only negligibly and have not considered it any longer in 
our calculations. The decay amplitude M. in (Al) then 
depends only on relativistic invariants Sij = (pi + pj) 2 , 
which are symmetric against transformation ip' 3 — > lix — 
ip' 3 . Using this symmetry when calculating the innermost 
integral in (Al) enables us to speed up the computing by 
a factor of two. 

In the diagrams with the a in the intermediate states, 
depicted in Figs. 1 and 2, also the product h ai<77T g cr7T7T 
plays a role (for definitions, see Sec. II B). We introduce 
additional free parameter 



which allows us to express that unknown product as 
a multiple of g ai aTrgt7-K-K, which is determined by the 
method described in Sec. II C. The parameter z itself 
will be obtained by minimalization of x 2 when fitting the 
experimental three-pion mass distributions in the three- 
pion decays of the r lepton. 

To be closer to experimental conditions, we do not 
calculate the (unnormalized) differential decay rate at 
certain values of the three-pion mass W or its square 
Q 2 = W 2 , but its averages over the experimentally given 
bins in W [7, 12] or Q 2 [9, 15]. 

When the contribution of the a[ to the a\ propagator 
is not considered, the calculated differential decay rates, 
and thus also the x 2 evaluated from them and data, de- 
pend on the following four parameters: (1) the nomi- 
nal a\ mass m ai , (2) the nominal a\ width r oi , (3) the 
aip-K Lagrangian mixing parameter sin0, and (4) the off- 
mass-shell coupling constant ratio z defined by Eq. (21). 
Quantity y (16) is not an extra parameter, condition (7) 
determines it as an implicit function of m ai and sin#. 

The ARGUS, OPAL, and CLEO experiments present 
the three-pion mass spectra in the acceptance corrected 
number of events. Both ALEPH spectra are normalized 
to the integrated branching fractions. As the outcome of 
our model is not normalized (the coupling of the a\ meson 
to the W boson is not fixed by meson dominance [64]), 
we opt to compare just shape of the mass distribution 
and introduce five multiplicative constants. The values of 
them are obtained by minimizing the individual x 2 's for 
each experiment while keeping the common parameters 
(m 0l , r oi , sin#, and z) fixed. 

To get a quick insight into the dependence of the 
quality of the fit on the Lagrangian mixing parameter 
sin0, we first fix the basic a\ parameters at the "stan- 
dard" values, frequently used in theoretical considera- 
tions, namely, m ai — 1.23 GeV/c 2 and T ai = 0.4 GeV. 




-1 -0.5 0.5 

sin 6 



FIG. 3. Dependence of \ 2 divided by the number of experi- 
mental points on the Lagrangian mixing parameter sin 9 for 
individual data sets. The interval of sin# suggested by the 
electron-positron annihilation into four pions [53, 54] is de- 
picted as a short abscissa in the bottom right corner. 



We also set z = and calculate the ratio of the usual 
X 2 to the number of experimental points N for each of 
the five data sets as a function of sin 9. The results are 
shown in Fig. 3. It is clear that the choice of the correct 
aip-K Lagrangian is of the utmost importance for obtain- 
ing a good agreement with the data. The fact that all 
five experiments point to the same narrow region in sin 9 
is extremely important. In addition, this region overlaps 
with the interval (0.457,0.471) based on the results of 
the model [53, 54] of the electron-positron annihilation 
into four pions built around the same Lagrangian (1). It 
indicates the soundness both of the present model and of 
the e + e~ annihilation model. 

The region sin 9 > 0.5 is not shown in Fig. 3 because for 
those values of sin# it is impossible to satisfy condition 
(7) by procedure described in Sec. II C. The square of 
parameter y, defined by Eq. (16), acquires negative, i.e. 
unphysical, values, which mean the negative branching 
ratio of decay (12). 

In the next step we allow all four parameters to vary 
and use the CERN computer library program MlNUlT of 
James and Roos [65] for finding their values that mini- 
mize x 2 f° r the three data sets defined in Sec. III. The 
results are summarized in Table II. As always, the as- 
sessment of errors of the parameters is a difficult task. 
We combined the errors provided by Minuit, which re- 
flect the errors of experimental data, with our estimates 
of the errors induced by the uncertainties of the input pa- 
rameters (the a mass, various coupling constants). The 
agreement of our model with the ALEPH 7r~7r°7r° data 
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TABLE II. Results of fitting various data sets. Only ai(1260) considered. Parameter y is defined by Eq. (16), parameter z by 
Eq. (21). All= ARGUS [7] + OPAL [9] + CLEO [12] + both ALEPH [15] data sets. For comparison, the values of Lagrangian 
mixing parameter from the e + e~ annihilation into four pions are also shown. 



Data 


Type 


x 2 / NDF 


C.L. (%) 


m ai (MeV) 


r oi (MeV) 


sin# 


V (GeV 2 ) 


z 


ALEPH [15] 




119.1/111 


28.25 


1220 ± 20 


418 ± 40 


0.460 ± 0.004 


0.094 ±0.010 


0.31 ±0.03 


ALEPH [15] 


- 
TV TV TV 


51.5/111 


100.00 


1256 ± 10 


443 ± 15 


0.466 ± 0.004 


0.111 ±0.022 


0.12 ±0.16 


All 


Mixed 


357.7/321 


7.74 


1232 ± 25 


431 ± 25 


0.463 ± 0.005 


0.099 ± 0.009 


0.30 ±0.05 




e + e~ -> 7r + 


7T~7T + 7r~ [53] 
7r~7T + 7T~ & 7T 


+7T-7TV 


[54] 




0.460 ± 0.003 
0.466 ± 0.005 







is perfect (confidence level 100%), the agreement with 
two other data sets is satisfactory by measures usually 
accepted in the high energy physics (x 2 /NDF The 
values of the a\ mass obtained from the three data sets 
are mutually compatible, as well as those of the ct\ width. 
What is especially remarkable are the values of the La- 
grangian mixing parameter sin 6. Not only their small 
errors (« 1%) and mutual consistence, but also a perfect 
agreement with the values obtained from the analyzes of 
the e + e~ annihilation into four pions. Parameter z char- 
acterizes the part of the aiair Lagrangian that acts only 
for virtual a\ and cannot be compared with anything yet 
(the e + e~ annihilation model [53, 54] did not consider 
the cr7r intermediate states). 

Now we add the a[ contribution (18) to the a\ propa- 
gator (5). Not to increase the number of free parameters 
too much, we fix the mass and width of the a[ at the 
PDG 2008 values 1647 MeV and 254 MeV, respectively. 
The same approach was used by the CLEO Collaboration 
[12], just their values were a little different (1700 MeV 
and 300 MeV). The number of the free parameters thus 
increases by two [the real and imaginary parts of a, see 
(18)]. The results of the x 2 minimalization procedure are 
shown in Table III for all three data sets. The compar- 
ison of Tables II and III shows that the addition of the 
a[ resonance to the a\ propagator greatly improves the 
agreement of the model with data in all cases. For the all- 
charged-pions ALEPH data [15], the x 2 drops from 119.1 
to 30.7 and the confidence level rockets from 28.25% to 
100%. The improvement of the confidence level is even 
more substantial for the third data set, where the to- 
tal x 2 is a sum °f the individual x 2 's for the ARGUS, 
OPAL, CLEO, ALEPH tt-jt+tt-, and ALEPH 7r-7r°7r° 
data. The mass and width of the a\ as well as other 
two free parameters (sin# and z) are very stable against 
the inclusion of a[. Their new values (Table III) differ 
only very little from the corresponding old ones (Table 
II). Also the values obtained from different data sets 
are mutually compatible. This is true also for two new 
parameters Re a and Im a. 

The calculated three-pion-mass distribution is com- 
pared to the 7r - 7r + 7r~ data of ALEPH Collaboration [15] 
in Fig. 4. We have just learned that the ai(1640) res- 
onance greatly improves the agreement with data. It 
is therefore a little surprising that there is no bump 
or shoulder corresponding to this resonance visible in 




0.5 1 1.5 2 2.5 

Q z (GeV) 



FIG. 4. Three-pion-mass-squared distribution calculated as- 
suming both ai and a[ contributions and compared to the 
ALEPH [15] t~ — > !/ T 7r~7r + 7r~ data. The model parameters 
taken from Table III. 



Fig. 4. To investigate this conundrum we calculate the 
model distribution in three cases: (1) both ai(1260) and 
ai(1640) terms in the a\ propagator (this is the curve 
presented already in Fig. 4); (2) only the ai(1260) term 
5 in ai the propagator; (3) only the a\ term (18). The 
model parameters in all three cases are identical. They 
are taken from the ALEPH ir~ir + ir~ row of Table III. 
The findings, see Fig. 5, show that the underlying mech- 
anism leading to the agreement with data is somewhat 
surprising. The final distribution is a result of the de- 
structive interference between the dominant amplitude 
containing the di(1260) propagator (5) and the ampli- 
tude containing the a[ propagator (18). Similar analysis 
performed for the t~ — > z/ T 7r~7r°7r° decay, see Fig. 6, 
leads to the same conclusion. 
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TABLE III. Results of fitting various data sets. Both ai(1260) and ai(1640) are considered. For definition of a, see Eq. (18). 

Data Type xV NDF C.L. m ai (MeV) r ai (MeV) sin<9 y (GeV 2 ) z Re a Im a 

ALEPH [15] 30.7/109 100% 1218 ± 19 418 ± 30 0.457(4) 0.106 ± 0.019 0.34 ± 0.03 -0.30 ±0.10 0.31 ± 0.06 

ALEPH [15] tt-ty°tt° 12.3/109 100% 1255 ± 18 455 ± 15 0.457(6) 0.148 ± 0.025 0.36 ± 0.14 -0.34 ± 0.13 0.29 ± 0.10 

All Mixed 219.5/318 100% 1233 ± 18 431 ± 20 0.459(4) 0.114 ±0.014 0.34 ± 0.05 -0.31 ±0.10 0.32 ± 0.09 




FIG. 5. Investigating the role of the a[ resonance in fitting 
the ALEPH tt~-k + -k~ data. Full curve: the complete calcu- 
lation shown in Fig. 4; Dotted curve: parameters unchanged, 
but only the <zi(1260) term (5) in the a\ propagator; Dashed 
curve: parameters unchanged, but only the a[ term (18) in 
the ai propagator. Note the change of scale against Fig. 4. 



FIG. 6. Investigating the role of the a[ resonance in fitting 
the ALEPH 7r~7r°7r data. Full curve (buried in data): the 
complete calculation; Dotted curve: parameters unchanged, 
but only the ai(1260) term (5) in the ai propagator; Dashed 
curve: parameters unchanged, but only the ai(1640) term 
(18) in the a\ propagator. 



V. SUMMARY AND CONCLUSIONS 

The main message of this study is that the form of 
the aip-K Lagrangian is the decisive factor for achieving 
a good model description of the three-pion decays of the 
tau lepton. This is again illustrated in Fig. 7 where the 
total x 2 j which is calculated as a sum of the individual 
% 2 's from the five experiments, is divided by the total 
number of experimental points and plotted as a function 
of the Lagrangian mixing parameter sin#. Two different 
cases are considered: (1) only ai(1260) included in the a\ 
propagator, (2) both ai(1260) and ai(1640) included. In 
contrast to Fig. 3, the other parameters are fixed at their 
optimal values taken from the appropriate tables (Tabs. 
II and III). Even if the curve (2) is shifted a little toward 
smaller values of sin#, the minima of both curves fall to 
the interval found in the model of the e + e~ annihilation 
into four pions [53, 54]. 

Our further finding, even not documented in this work 
in detail, concerns the form of the a\ propagator. We 



have found that the running mass form (5), suggested 
and already used in several papers [12, 18, 19], provides 
a better fit to the taon three-pion decay data than simpler 
forms with a constant a\ mass and a constant or energy 
dependent a\ total decay width. 

The a\ running mass squared is given by the dispersion 
relation. In this work we have chosen a once-subtracted 
version (9). As the input for the dispersion relation, the 
energy dependent total decay width of the a\ for all s 
above the three-pion threshold is required (9). We ap- 
proximated it as a sum of the decay widths to the three 
pion final states (via the irp and ttu intermediate states) 
and the KKn final states (via K*(892)K + ex.). A typi- 
cal behavior of the energy dependent total width is shown 
in Fig. 8. The hump centered around i/s w 1 GeV devel- 
ops as the mass of the two-pion subsystem falls predom- 
inantly first on the ascending and then on the descend- 
ing side of the rho propagator. We ignored the channels 
P(1450)tt, / (1370)tt, and / 2 (1270)tt, which have been 
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FIG. 7. sin 8 dependence of the sum of x s from all five exper- 
iments divided by the total number of experimental points. 
The other parameters are kept at the optimal values from the 
"All" row of the corresponding tables. Dotted curve: only 
ai(1260) (Table II) included in the oi propagator; Full curve: 
both ai(1260) and ai(1640) included (Table III). The range 
of sin 6 from the electron-positron annihilation into four pions 
[53, 54] is shown as a short abscissa. 



0.8 




FIG. 8. Energy dependent width of the ai(1260) as a function 
of yfs. Parameters m ai , T ai , and sin# taken from Table III, 
row "All". 



a 1.26 









m a = 1.233 GeV 
T a = 0.431 GeV 
sin = 0.459 




, , , , 



,1 



LIB 

0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 
Vs (GeV) 

FIG. 9. Running mass of the ai(1260) as a function of ^/s for 
the same parameters as Fig. 8. The cross marks the point in 
which condition (6) is satisfied. 



seen in the a\ decays [39] and which open at higher s. 3 
This is probably the reason why the running mass be- 
haves wildly, see Fig. 9, and does not have a nice plateau 
around the nominal mass, as it did in the case of the 
p(770) [59]. 

Another important ingredient of our model are the na 
intermediate states. On one side, they enter the calcula- 
tion of the total decay width of the a\ resonance, which 
is necessary for constructing the running mass propaga- 
tor (5). On the other side, they contribute to the decay 
rates of the three-pion decays of the tau lcpton, Figs. 1 
and 2. 

To investigate the role of the ira intermediate states in 
the evaluation of the r~ — > z/ T 7r~7r°7r° decay width, we 
split the distribution depicted in Fig. 6 by the full curve 
into its irp and na components. The result is shown in 
Fig. 10. It is obvious that the it<j intermediate states 
play a unique role in describing the behavior of the dif- 
ferential decay width at small three-pion masses. What 
is a little suspicious, is the large magnitude at the inter- 
mediate masses. To see whether it is reasonable or not, 
we integrate the distributions to get the branching ratio 



B 



r(. 



z/ T 7r-7r°7r°) 



with the result B « 41%. This number is more than twice 
higher than the experimental value of (16.18 ± 3.85 ± 



3 The inclusion of them would bring additional free parameters, 
what we wanted to avoid. 
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FIG. 10. Decomposition of the three-pion mass squared dis- 
tribution into the contributions from the up and no interme- 
diate states. The parameters were taken from the ALEPH 
Tr-Tr'V row of Table III. 



ing parameter sin 9 perfectly agrees with the value ob- 
tained from the e + e~ annihilation into four pions. 

(2) Our confirmation of the existence of the di(1640) res- 
onance with the mass and width compatible with the 
PDG [39] values is based on the increase of the confi- 
dence level from 7.7% to 100% after the ai(1640) has 
been included. 

(3) We have explained why the a\ (1640) resonance, which 
is important for getting a good agreement with data, is 
not visible in the three-pion-mass spectrum as a bump 
or shoulder. 

(4) From the common fit to the data from five experi- 
ments we have obtained the following results: 

Mass of the en (1260) m 0l = (1233 ± 18) MeV; 
Width of the ai(1260) T Ql = (431 ± 20) MeV; 
Lagrangian mixing parameter sin 6* = 0.459 ± 0.004. 
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1.28)% obtained by the CLEO Collaboration (Rcf. [12], 
Table III). 

The source of this obvious deficiency of our model is 
the following: The important parameter y, which reg- 
ulates the rate of the a\ — > ttu — > 3tt transition, was 
obtained from the condition that the derivative of the 
running mass at the nominal-mass point should vanish 
(7). But the absence of higher decay channels, which 
influence the values of the running mass at all s, may 
modify the resulting value of y significantly. The larger 
than correct y may mimic the missing channels. 

A new feature of our work, which, to our knowledge, 
has not appeared in the literature yet, is that we fit the 
data from several experiments simultaneously. We in- 
tend to continue in this approach and include not only 
the data concerning the three-pion decays of the tau lep- 
ton, but also the experimental results from other weak, 
electromagnetic, and perhaps also strong interaction pro- 
cesses. The natural candidate is the electron-positron 
annihilation into four pions, for which a model based on 
the same Lagrangian as here has already been built. The 
value of the Lagrangian mixing parameter we have ob- 
tained here perfectly agrees with values obtained from 
the e + e~ annihilation into four charged pions [53] and 
from the combined fit to both annihilation channels [54] . 

To summarize: 
(1) We have shown that the right form of the aipn La- 
grangian is extremely important for obtaining a good 
agreement with data. We have obtained an unprece- 
dented confidence level of 100% for all three sets of data 
we considered. The optimal value of the Lagrangian mix- 



Appendix A: Differential decay rate formula 

We use the following formula for the differential de- 
ca y rate in the invariant three-particle mass W = 
y/ (P2 + P3+ Pi) 2 in a four-body decay a — > 1 + 2 + 3 + 4: 

dT rW-m 2 

/l p\ /"27T 

rfcos6>2 / dcos6' 3 d</?3 |A4| 2 . (Al) 
-1 J-i Jo 

The asterisk denotes the (2,3,4) rest frame, the prime the 
(3,4) rest frame. 771,34 is the mass of the system consisting 
of particles 3 and 4, E 34 = E^+E% and Pg 4 = vl + Vl = 
— P2 are its energy and momentum, respectively, in the 
(2,3,4) rest frame. In the rest frame of the parent particle 
a the momentum of particle 1 points along the negative z- 
axis. In the (2,3,4) rest frame, the momentum of particle 
2 lies in the (xz) plane. 

Appendix B: Integrating over a narrow peak 

Let us assume that we need to evaluate an integral over 
an interval that includes a narrow resonance peak 

f S2 f( s ) 
Q = J S1 ( S -M2( s ))2 +TO 2 r2(s) ds > ( B1 ) 

where f(s) is a slowly varying function. Further, let 
the two functions in the denominator satisfy conditions 
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M 2 (m 2 ) — m 2 and T(m 2 ) = 7. If 7 <C m then the in- 
tegrand is rapidly varying function of s and a numerical 
quadrature of very high order is required to get reliable 
results. After introducing a new variable £ by substitu- 
tion s = m 2 + TO7tan (c £ + d), where c = (a 2 — ai)/2, 
= (ai + d2)/2, ai = arctan{(si — m 2 )/(m7)}, and 
a 2 = arctan{(s2 — to 2 )/(to7)}, the integral (Bl) becomes 



which can be safely evaluated using, e.g., the Gauss- 
Legendre quadrature. We apply this method for calcu- 
lating the integrals containing the square of the if* (892) 
propagator (20). In that case M(s) = m = vtik* ■ 



Q 



a 2 - a 1 



(s — m 2 ) 2 + m 2 7 2 



2m 7 y_! (s - M 2 (s)) 2 + m 2 r 2 (s) 



/(*) , 
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